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0\ ' Abstract 






The standard General Relativity results for precession of particle orbits and 
for bending of null rays are derived as special cases of perturbation of a 
^^ I quantity that is conserved in Newtonian physics, the Runge-Lenz vector. First 

this method is applied to give a derivation of these General Relativity effects 
for the case of the spherically symmetric Schwarzschild geometry. Then the 
pg ■ lowest order correction due to an angular momentum of the central body is 

OO . considered. The results obtained are well known, but the method used is 

rather more efficient than that found in the standard texts, and it provides a 
good occasion to use the Runge-Lenz vector beyond its standard applications 
in Newtonian physics. 
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I. INTRODUCTION 

Light bending and perihelion precession are the two most important effects on orbits 
caused by the General Relativity corrections to the Newtonian gravitational field of the sun. 
The standard derivation treats these two effects in different ways, without any apparent 
connection between them. Yet, in the usual Schwarzschild coordinates they are both due to 
the same, single relativistic correction to the Newtonian potential, so it is of some interest 
to use the same method to derive either effect. 

The key to the present unified treatment is the Runge-Lenz vector. In Newtonian physics, 
where the two effects are absent, this vector is constant and points from the center of 
attraction to the orbit's perihelion.^] Its non-constancy in General Relativity therefore is a 
measure of either effect. The Runge-Lenz vector was established as a useful tool by 1924 at 
the latest, but it did not become popular until the 196O's.0 

Since then a number of papers that exploit its advantages have graced the pages of this 
Journal (see ref. and the references cited therein), and the results to be reported here can 
in essence be found in earlier papers, but the unified viewpoint vis a vis General Relativity 
is perhaps new. In addition the "magnetic" gravitational effects due to a rotating central 
body are treated here with this method. 

II. GENERAL RELATIVISTIC EQUATIONS OF MOTION (NO ROTATION) 

The motion to be considered is that of a "test particle" that moves along a geodesic in the 
spacetime exterior to the central body. If this body is non-rotating, spherically symmetric, 
and has total mass M, the exterior spacetime geometry is described by the Schwarzschild 
line element 

ds^ = -[l \dt^ + :^ + r^dn^. (2.1) 

r 

Here the coordinates t, r, 9 and cf) are one of many equally valid choices for labeling spacetime 
points, but they can nevertheless be invariantly characterized [§]. (For example, 47rr^ is 
the area of the sphere r = const, t = const, and d/dt is a timelike Killing vector.) The 
geodesic law of motion is essentially equivalent to the conservation laws that follow from the 
symmetries of the geometry and the conservation of rest mass. Because of the principle of 
equivalence we may assume without loss of generality that the rest mass is e, where e = 1 for 
particles of finite rest mass, and e = for light (photons). The other conserved quantities 
correspond to angular momentum L, and energy E. We use r to denote the proper time 



^Although our treatment is not confined to the solar system, we use this term to denote the 
point of closest approach to the center because it seems more familiar (and more etymologically 
consistent) than the more correct term, pericenter. 

^For a history of the Laplace- Runge-Lenz vector, see ref. ffl. 



along the geodesic.^] As usual we can choose the vector L to be normal to the plane 9 = 7i/2 
and we then have^ 

.^.f P.3) 

Equation (2.4) has the form of conservation of energy in an effective potential. Except for 
the presence of r instead of t, Eqs (2.2 - 2.4) are the same as the Newtonian equations of 
motion of a particle of unit mass and total energy f^ in a potential V = —eM/r — ML'^/r^. 
Thus for particles as well as for light the relativistic motion in proper time r is the same as 
the Newtonian motion in Newtonian time t if the potential is modified by the single term 
—ML'^/r^. We note that no slow motion assumption or other approximation is involved 
in this correspondence. If we are only interested in the orbit equation, then the difference 
between t and r does not matter, because either one will be eliminated in the same way in 
favor of via Eq (2.3). 

III. SECULAR CHANGE OF ORBITS 

We treat the modification ML'^/r^ of the Newtonian potential as a perturbation^ and 
compute the consequent changes in direction of the Runge-Lenz vector, defined by 

A = V X L - eMer (3.1) 

where all bold-face quantities are 3-vectors, and e^ denotes a unit vector in the r-direction. 
The time parameter is r as above, so that, for example, v = dr/dr. The rate of change of 
A is (as, for example, in ref. p) 



r'"l^-^-Mr^^r^r^e,. (3.2) 



rfA _ / ^dV_ _ ^\de^_ (SMI^ 

dr \ dr J dr \ r"^ ) dr 

The direction of A therefore changes with angular velocity 



A X A /3ML2\ d(h ^ , , 



^For the case of light, r is an affine parameter that is defined only up to scale transformations. The 
quantities E and L are therefore similarly defined only up to such re-scaling. The final, physical 
results will contain only ratios of such quantities. 



^For a derivation, see references |^,^ or section IV below. 

^This means that M/r <C 1, where r is a typical orbit radius, which follows if we assume M/L ^ 1 



and its total change when the particle moves from 0i to 02 is (assuming this change is small, 
and that A is originally in the = direction) 

Aa=f''u;dr = 3ML^f'''-^. (3.4) 

When A is constant, and in the direction = 0, we have 

A • r = Ar cos = L^ - eMr. (3.5) 

The bound orbits (e = 1) are therefore ellipses with eccentricity e = A/M and semi-major 
axis a = j^,^_^2\ - For the unbound, straight orbits of light rays (e = 0), L'^ / A = 6 is the 
impact parameter, and because these orbits are traversed at the speed of light we have 
L/ E = b. In either case A points from the center to the perihelion. When A changes slowly 
the orbits still have this approximate shape, but their orientation and shape will change 
slowly. We calculate the lowest order changes due to the General Relativistic correction in 
Ves by substituting the unperturbed orbit (3.5) into Eq (3.4): 

3M /•<^2 



AL^ 



r<p2 
/ (Acos0 + eM)2cos0rf0. (3.6) 



A. Perihelion precession 

For a particle in a bound orbit it is customary to find the angular change of the perihelion 
during one revolution (in 0) of the particle. Because A points to the perihelion, this angle 
is given by Eq (3.6), when 02 — 0i = 27r. 

3M /•2'r (Acos0 + M)2 , ,, GttM^ QttM , „, 

^^ = uk — -A — ''''''^^^=^r = w^y ^^^ 

This is the usual perihelion formula. 

B. Light bending 

Here also the deflection is given by Aa of Eq (3.6), but changes from — 7r/2 to 7r/2 
with respect to the perihelion (and of course e = 0): 

3M rVa 4MA 4M , , 

This is the usual light deflection formula. 

That the light deflection should follow from the same, 0(l/r^) correction to the New- 
tonian effective potential as the perihelion rotation may be somewhat surprising, because 
the deflection is frequently heuristically explained as an action of the Newtonian 0(l/r) 
potential on light. Indeed, an effective potential for dr jdt would contain 0(l/r) terms, but 
in the present choice of variables these are absent — illustrating once again the arbitrary 
nature of coordinates in a generally covariant theory. 



IV. SLOWLY ROTATING CENTRAL BODIES 

If the body is slowly rotating^ in the 0-direction with angular momentum J, the metric 
(2.1) is modified by the Lense-Thirring term [Q, —— siv? Odcpdt and by a quadrupole term 
that describes the distortion of the body. The quadrupole term can be treated in the 
same way as the relativistic correction, and will therefore not be further considered. The 
Lense-Thirring term breaks the spherical symmetry, so on symmetry grounds only p^, the 
angular momentum conjugate to (called L below), is conserved. Nonetheless the "total 
angular momentum" Q"^ = Pg + cot^ 6'p^, where pe is the momentum conjugate to 6', is also 
conserved to first order in J . Thus the entire motion can be formulated in terms of conserved 
quantities, and one finds that the Newtonian angular momentum L precesses around the 
2;-direction.[] However, we will confine attention to the case when J and L are parallel, 
the motion is confined to the equatorial plane, and the general relativistic correction is 
completely described by the behavior of A. 

A. Equations of Motion and Runge-Lenz vector 

A Lagrangian C for a particle or a light beam moving the in equatorial plane of metric 






ds'^ = ~ (l - —] dt^ + -^ 2M + ^^^^^ - — sin^ Od(t)dt (4.1) 



is given by 



2M\ .o r^ ,., 4J 

1 2M 



2C = -[l- — ]i^ + - 2M+^V'-— <^^- (4-2) 



Here the dot indicates differentiation with respect to proper time; this implies that C is 
constant, 

2C = -e. (4.3) 

Since C is independent of t and we have the conserved quantities, 

^ dC ( 2M\ . 2J . . dC 0-, 2J. ,_, 

-E= —^ = -{1 ]t L = ^ = r^(p 1. 4.4 

dt \ r J r ^ d(p r ^^ 



^We work only to first order in J; more precisely, we assume J ~ ML and, as before, M/r ^^ A'P/L^ ~ 
£ < 1, so that J/r'^ -- e^/^. 

The integration of the equations of motion for orbits of general orientation was one of the aims 
of the Lense and Thirring paper Q. For a treatment using the Runge-Lenz vector, see references 
|0]. For a summary of all the relativistic effects on orbits see reference [§. For the geometrical 
reason for the conservation of Q^ see reference S . 



We solve for t and (f) to first order in J, 

• L 2JE . E 2JL ,, ^, 

r 

and substitute into Eq ( [4.3|) to obtain a "conservation of energy" in an effective potential, 



. , 2eM L^ 2ML2 4JLE , , 

E^-e = f^ + — — + — . 4.6 



The effective potential in this equation contains two non-Newtonian terms. The first was 
already encountered in Eq ( |2.4| ) and causes the "standard" relativistic correction; the second 
is due to the Lense-Thirring addition to the metric.^ 

Because in the rotating case there is a difference between kinematic angular momentum 
(r^0) and canonical angular momentum L, the Runge-Lenz vector A can be defined in 
various ways, but there is no difference in the precession rate one calculates from them. A 
convenient choice is 

A = V X ("l - ^] - eMer = (y - ^m\ e, - r (l - ^) e^ (4.7) 

because it simplifies the equation of motion for A, and still gives elliptical orbits for any J 
when A is constant, 

A-er = Acos0 = eM. (4.8) 

r 



The equation of motion for A can be derived from Eqs ([4.5| 

. /3ML2 SeMJE 2JE(e-E^)\- , ^ ^^ 

By substituting Eq (|4.9|) into Eq ( p.3|) and integrating, using (^l8| ) for 1/r, we now find that 
the total change in A when the particle moves from 01 to 02 is 

J (tp^^^°^'^ + ^^) — ^(Acos0 + eM) + ^— ^Jcos0c?0. 

(4.10) 



^With our assumptions as spelled out in footnote 6 the Newtonian terms of the effective potential 
are of order e, both non-Newtonian terms are of order e^, and typical terms that are neglected are 



B. Perihelion motion 



To obtain the perihelion motion we evaluate Equation ( |4.1CI| ) over one revolution 



^1 



0, 02 = 27r) with e = 1. Since the particle-velocity is nonrelativistic, we may set ii^ = 1 to 
the lowest order :0 

_ GttM^ SnJME _ GttM 8nJ 

""~^^ I^~ ~ a(l - e2) ~ MV2 (a(l - e^)f' ' ^ '^^^ 

The first term is the "standard" general relativistic precession already found in Section III, 
and the second term is due to the rotation of the central body.Q For nearly circular orbits 
we can interpret this second term as due to two causes: one is the rotation in (p of the 
"locally non-rotating observer" that makes the Lense-Thirring term of Eq ( f4.2| ) disappear 
at the radius of the particle, an amount ATcJ/aL; the other is the "differential rotation" due 
to the 1/r^ fall-off of the second non- Newtonian term in the effective potential of Eq ( [4.6p . 
This contribution causes precession by an amount — 127rJ/aL, in the same way as the first 
non-Newtonian term causes the "standard" precession.[^ 

C. Light bending 

For the effect on light we put e = in Eq ( |4.1(]| ) and integrate from = — 7r/2 to 7r/2 as 
in section |T| B, 

3MA 4 2JE^ AM AJ , , 

The effect of the angular momentum J on both precession and bending is negative .0 This 
is the same "differential" dragging effect that makes a gyroscope in the equatorial plane 
precess in the opposite direction to the central body's rotation [Q. 

V. CONCLUSIONS 

For non-rotating spherically symmetric central masses we have seen that the two impor- 
tant general relativistic corrections to the Newtonian gravitational motion, namely perihelion 



^Formally, this follows from the M^ <^ L^ assumption and requiring bound orbits. 

^'^The last term in Eq ( [4.1l| ) changes sign if L is antiparallel to J. Both terms are of order e. 

^^For non-equatorial orbits the first contribution is a precession about J, whereas the second 
contribution is a precession about L, proportional to J • L. 

^^However, the relative the contribution of J to light bending is less: we have M/b ~ e, but 



precession and light bending, follow from the same correction term in the effective radial 
potential; and that either effect can be viewed as change in the Runge-Lenz vector associ- 
ated with the orbit. Because both effects follow by simple evaluation from one formula (Eqs 



4.10|) , the effort is only about half of the usual procedure; moreover it gives occasion to 
review and apply the Runge-Lenz vector. We have shown the extension of this calculation 
to equatorial orbits of a rotating body; relativistic corrections to parabolic and hyperbolic 
orbits can similarly be evaluated by this method. 
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